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ABSTRACT 

In the presence of a magnetic field, rotational energy can be extracted from black holes via the 
Blandford-Znajek mechanism. We use self-similar advection dominated accretion (ADAF) models to 
estimate the efficiency of this mechanism for black holes accreting from geometrically thick disks, in the 
light of recent magnetohydrodynamic disk simulations, and show that the power from electromagnetic 
energy extraction exceeds the accretion luminosity for ADAFs at sufficiently low accretion rates. We 
consider the detectability of isolated stellar mass black holes accreting from the ISM, and show that 
for any rapidly rotating holes the efficiency of energy extraction, eez = L-qz/Mc 2 , could reach 10 -2 . 
The estimated total luminosity would be consistent with the tentative identification of some EGRET 
sources as accreting isolated black holes, if that energy is radiated primarily as gamma rays. We discuss 
the importance of emission from the Blandford-Znajek mechanism for the spectra of other advection 
dominated accretion flows, especially those in low luminosity galactic nuclei. 

Subject headings: accretion, accretion disks — black hole physics — magnetic fields — galaxies: nuclei 
— gamma rays: observations — X-rays: galaxies 



1. INTRODUCTION 

The formation of both stellar mass and supermassive 
black holes could well lead to rapidly rotating Kerr black 
holes. If so, this reservoir of rotational energy constitutes 
a potentially important energy source for driving jets and 
emission from black hole systems. Moreover, this energy 
can be extracted if the black hole is embedded in a mag- 
netic field (Blandford & Znajek 1977). Since all black 
holes of current astrophysical interest are probably ac- 
creting from magnetized disks, this has led to suggestions 
that the Blandford-Znajek process plays a vital role in 
Active Galactic Nuclei (AGN) and other accreting black 
hole systems (Begelman, Blandford & Rees 1984; Wilson 
& Colbert 1995; Moderski, Sikora & Lasota 1998; Paczyn- 
ski 1998; Lee, Wijers & Brown 1999). 

Recent papers (Ghosh & Abramowicz 1997; Livio, 
Ogilvie, & Pringle 1999; Li 1999) have re-evaluated the 
importance of the Blandford-Znajek mechanism for AGN 
in light of advances in the understanding of accretion disks 
- in particular, the demonstration that turbulence driven 
by magnetohydrodynamic instabilities (Balbus & Hawley 
1991) leads to disk angular momentum transport (Hawley, 
Gammie & Balbus 1995; Brandenburg et al. 1995; Stone 
et al. 1996; for a review see e.g. Hawley & Balbus 1999). 
Two principal effects lead to substantial downward revi- 
sions in the estimated efficiency of the Blandford-Znajek 
mechanism for thin disks. First, in an accretion disk both 
the viscosity and the magnetic diffusivity are generated 
by the same turbulent processes, and so have comparable 
magnitudes. This implies that any external magnetic field 
can diffuse outward through the inflowing gas, and will not 
accumulate on the hole. Second, the strongest magnetic 
fields generated by a thin disk are expected to have char- 
acteristic scales of order the disk scale height H, with only 
weaker large scale fields threading the black hole and con- 
tributing to the extraction of the rotational energy. Both 



arguments are weakened for thick disks for which H ~ R. 

In this Letter, we consider the Blandford-Znajek contri- 
bution to the luminosity of black hole systems at extremely 
low accretion rates. Motivated by the hard high energy 
spectra and low luminosity of such systems, Narayan & 
Yi (1994; 1995) have developed an advection dominated 
accretion flow (ADAF) model, in which radiatively ineffi- 
cient accretion occurs in a hot, geometrically thick geome- 
try (see also Rees et al. 1982; Ichimaru 1977). Strong disk 
winds may also contribute to the observed dearth of emis- 
sion (Blandford & Begelman 1999). As Rees et al. (1982) 
noted, the low radiative efficiency cadaf = L adaf / Mc 2 
of ADAFs means that even rather weak emission arising 
from the Blandford-Znajek mechanism could make an im- 
portant contribution to the overall luminosity. Here, we 
concentrate on the extreme case of isolated black holes ac- 
creting from the interstellar medium (ISM), and show that 
the power from the Blandford-Znajek process substantially 
improves the prospects for detecting any rapidly rotating 
holes. We also comment on the implications for other very 
low luminosity accreting black holes, particularly those in 
quiescent galactic nuclei. 

2. EVALUATING THE BLANDFORD-ZNAJEK POWER 

For a black hole of mass M and angular momentum J, 
with a magnetic field B± normal to the horizon at Rh, 
the power arising from the Blandford-Znajek mechanism 
is given by (e.g. Ghosh & Abramowicz 1997; Macdonald 
& Thorne 1982; Thorne, Price & Macdonald 1986), 
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where J max = GM 2 /c is the maximal angular momentum 
of the hole. The factor oj f = ilp(ilh — ftp)/^ depends on 
the angular velocity of field lines f2 f relative to that of the 
hole, fi/j. Conventionally, we will assume that top = 1/2, 



2 
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which maximizes the power output (Macdonald & Thorne 
1982; Phinney 1983; Thorne, Price & Macdonald 1986). 

To estimate Lbz we need to relate B± to the conditions 
at the inner edge of an ADAF. First though, we briefly 
consider the possibility that B± at the black hole hori- 
zon could exceed the disk field at the inner disk edge, if 
the disk is able to advect an organized external magnetic 
field inwards. Whether this is possible depends on the 
relative timescales for inward advection, controlled by the 
disk shear viscosity v, versus outward diffusion, controlled 
by the magnetic diffusivity rj. For a disk with magnetic 
Prantl number V = v/r), the ability of the disk to advect 
field lines inwards depends on the dimensionless param- 
eter T> = (R/H)V, where H is the disk scale height at 
radius R. Lubow, Papaloizou & Pringle (1994) found that 
"D <^ 1 was a necessary condition for strong inward field 
dragging (see also Reyes- Ruiz & Stepinski 1996). Since it 
is generally believed that in a turbulent medium V — 1 
(e.g. Parker 1979), this result implies that for a thin disk 
T> ~ (R/H) ^> 1, and external fields cannot accumulate at 
small radii. For thick disks, where H ~ R, the situation 
is less clear, but for a lower limit on B± we will assume 
that even in thick disks, advection of the field is ineffec- 
tive. The field at the black hole then reflects only the local 
disk magnetic field strength. 

To estimate the physical conditions in the inner part 
of the disk, we employ the vertically averaged, self-similar 
ADAF solution of Narayan & Yi (1994; Spruit et al. 1987). 
This solution assumes a Shakura-Sunyaev (1973) prescrip- 
tion for the shear viscosity, v — ac s H , where a is a free 
parameter and c s is the local disk sound speed. The char- 
acter of the solution depends on the parameter e', where 
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and 7 is the ratio of specific heats. The factor / mea- 
sures the efficiency of radiative cooling, so that / = 1 
corresponds to the limit of no cooling at all, while / = 
resembles the case of a thin disk in which cooling is effi- 
cient. For a <C 1, as suggested by numerical simulations, 
the solution for the sound speed and the radial velocity vr 
is (Narayan & Yi 1994), 
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where v% = {GM / R) 1 / 2 is the Keplerian velocity. Making 
use of H / R — c s /vk, and the continuity equation, 

p = -M/A-kRHvr, (4) 
we find that the pressure P = pc 2 is given by, 
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Studies that have dropped the requirements of self- 
similarity and vertical averaging have found that this 
solution provides a good approximation to ADAF flows 
(Narayan & Yi 1995). 

Numerical simulations of thin, magnetized accretion 
disks (Hawley, Gammie & Balbus 1995; Brandenburg et 



al. 1995; Stone et al. 1996), have shown that angular 
momentum transport is dominated by Maxwell stresses. 
The magnetic contribution to a typically exceeds the fluid 
stresses by an order of magnitude, so that a ~ Q! ma gnotic- 
The dominant field component is toroidal, with satura- 
tion occurring when P mag <C P. Initial results from global 
simulations (Armitage 1998) suggest that Narayan & Yi's 
assumption - that in thick disks the magnetic fields pos- 
sess most of their power in long wavelength modes of scale 
~ R - is reasonable. 

We assume that P mag = B 2 /8n ~ aP in the inner disk, 
and make the important assumption that B± w B. We 
further assume that the energy extracted from the hole is 
not all dumped into the internal energy of the accreting 
gas and lost invisibly across the horizon - the fate of the 
bulk of the accretion energy (Narayan, Garcia & McClin- 
tock 1997). This is reasonable provided that some fraction 
of the field lines threading the hole are open, for example 
as a result of a wind or a jet. Evaluating the Blandford- 
Znajek power for the case of a maximally rotating hole, 
where J = J ma x and Rh = GM/c 2 , we then obtain, 
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This corresponds to a constant efficiency, eez = Lbz/Mc 2 , 
given for the range of thick disk models by, 



14 



192 

VTo 

192~ 



for e = 1 
for e' = 
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This efficiency drops rapidly for more slowly spinning 
holes, due to both the (J/J m ax) 2 factor in equation (|j), 
and because the pressure in the disk, which should be 
evaluated at the marginally stable orbit, drops as R~ 5 ^ 2 . 
Nonetheless, it suggests that the Blandford-Znajek mech- 
anism will be an important power source in the overall 
energy budget of ADAFs around rapidly rotating black 
holes, with efficiencies of up to ~ 10 -2 . We note that if 
e' > 1 (corresponding to a thin disk), these scalings suggest 
still higher efficiencies. More realistically, the suppressing 
effects identified by Livio, Ogilvie & Pringle (1999) are ex- 
pected to be important in this regime. Of course, if a thin 
disk existed the holes would be readily detectable due to 
the high radiative efficiency (~ 0.1) of thin disk accretion. 
Additional electromagnetic extraction of energy from the 
inner disk, which Livio, Ogilvie & Pringle (1999) argue is 
likely to exceed that from the hole itself, would also boost 
the efficiency. 

3. ISOLATED BLACK HOLES 

The known neutron stars and stellar mass black holes 
in mass transfer binaries must be greatly outnumbered by 
the isolated populations of such remnants, accreting only 
at a very low rate from the ISM (e.g. Blaes & Madau 
1993). To date, only a few isolated neutron stars have 
been identified (Walter & Matthews 1997; Neuhauser & 
Trumper 1999). Identification of the corresponding black 
hole population would provide both a record of massive 
star formation in the galaxy, and insight into the behavior 
of accretion flows at extremely low M. 
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The accretion rate onto a black hole accreting from the 
ISM is given by the Bondi-Hoyle rate (Hoyle & Lyttleton 
1939; Bondi & Hoyle 1944). For a hole moving with veloc- 
ity Voo through a uniform medium with sound speed Coo, 
the accretion radius is, 
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Gas falling within a cylinder of this radius is accreted, so 
that the accretion rate is Msoudi = ^RaPoo^/voo + c 
where poo is the density far upstream of the accretor 

The most luminous black holes are expected to be those 
accreting from nearby, dense, molecular clouds (Fujita 
et al. 1998). For gas at temperatures of T ~ 10 2 K, 
Voo Coo, and the accretion rate is, 
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where we have taken a density Uoo (in particles / cm 3 ) 
typical for a molecular cloud. The accretion rate will be 
much lower for holes accreting from hotter phases of the 
ISM. For gas at Uoo ~ 1 cm -3 and T ~ 10 4 K, the same 
hole accreting transonically (voo w w 15 kms" 1 ) will 
have an accretion rate, Msondi = 8x 10 12 gs -1 . 

Some black holes in mass transfer binaries are believed 
to be rotating rapidly (Zhang, Cui & Chen 1997), and we 
assume that some fraction of isolated black holes also have 
rapid (J ~ J max ) rotation. Taking the Blandford-Znajek 
efficiency calculated earlier for the e' = 1 ADAF model 
(appropriate for a fully advection dominated, 7 = 4/3 
gas), the estimated power is, 
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1.3 x 10 35 ergs" 1 



("oo = 10 2 cm 3 , v. 



= 10 kms" 1 ) 
L BZ ^ 1-4 x 10 32 ergs" 1 
(noo = 1 cm" 3 , Woo = Coo = 15 kms" 1 ) 



(11) 



These luminosities are comparable to those predicted on 
the basis of models of spherical accretion of magnetized gas 
(Heckler & Kolb 1996; Ipser & Price 1982, 1983). Spheri- 
cal models are unrealistic at least for black holes accre ting 
at high Mach numbers. The luminosities in equation (11) 
are considerably higher than estimates based on the radia- 
tive efficiency of an ADAF (Fujita et al. 1998), which is 
very low for small M (Narayan & Yi 1995; Mahadevan & 
Quataert 1997). We note that at these accretion rates, the 
spindown time for the hole, even including the Blandford- 
Znajek losses, remains much greater than a Hubble time 
(e.g. King & Kolb 1999). 

Our estimates suggest that the prospects for detecting 
isolated black holes are substantially better than would 
be inferred on the basis of the low radiative efficiency of 
an ADAF, if a reasonable fraction of the holes are rapidly 
rotating. Indeed, these estimates are energetically con- 
sistent with associating isolated black holes with some 
unidentified EGRET gamma ray sources (Lamb & Ma- 
comb 1997; Dermer 1997; Yadigaroglu & Romani 1997). 



Dermer (1997) estimates that these sources have charac- 
teristic luminosities, in the 100 MeV - 5 GeV band, of ~ 
2.5 x 10 35 ergs -1 (for low latitude sources with |6| < 10°), 
and — 6 x 10 32 ergs -1 (for sources with |6| > 10°). This 
would be roughly consistent with the estimates for rotating 
lOM© black holes, provided that the efficiency of gamma 
ray emission was very high. High energy emission in AGN 
occurs via similar processes (e.g. Levinson & Blandford 
1996). We note that the rough proportionality of the lu- 
minosity of the unidentified sources on the ISM density 
(Dermer 1997) would be consistent with the approximately 
constant efficiency of the Blandford-Znajek mechanism. It 
would not be consistent with an ADAF model, whose ra- 
diative efficiency declines at low M (Narayan & Yi 1995). 

4. ADVECTION DOMINATED FLOWS IN GALACTIC 
NUCLEI 

Similar considerations apply to other low accretion rate 
flows that are advection dominated. The radiative effi- 
ciency of an ADAF is given approximately by, 



eadaf = 0.1 



(12) 



where m = Af/A/g^, and the Eddington accretion rate for 
an electron scattering opacity k cs , and fiducial efficiency 
Ves = 0.1, is defined as A/Edd = ^GM/r)ESK es c (Narayan 
& Yi 1995). For an estimated a = 0.1, this implies that 
the Blandford-Znajek effect could be an important con- 
tributor to the overall power output for rh <J 2 x 10 -3 . 

In addition to isolated black holes in the galactic disk, 
some supermassive black holes in quiescent galactic nuclei 
also fall into this regime (Fabian & Rees 1995). The weak 
emission from the nuclei of these galaxies provides both a 
test of the ADAF model (Di Matteo et al. 1999a), and 
a possible contributor to the X-ray background (Di Mat- 
teo & Fabian 1997; Di Matteo et al. 1999b). Taking our 
estimates at face value, the good agreement of computed 
ADAF spectra with observations in some systems - espe- 
cially the galactic center (Mahadevan 1998) - would be 
consistent with the black holes in those systems not being 
close to maximally rotating. 

5. DISCUSSION 

In this Letter, we have discussed the importance of 
the Blandford-Znajek (1977) mechanism for ADAFs at 
low accretion rates. For these thick disks, which are 
similar to the tori originally considered by Rees et al. 
(1982), the arguments advanced for the likely irrelevance 
of the Blandford-Znajek mechanism (Ghosh & Abramow- 
icz 1997; Livio, Ogilvie & Pringle 1999) are considerably 
weaker than for thin disks. Even with pessimistic assump- 
tions, the total electromagnetic power liberated from the 
hole or inner disk will exceed the radiative luminosity for 
sufficiently low accretion rates. For maximally rotating 
black holes, and self-similar ADAF models, we obtain the 
criteria m = M/M^dd < 2 x 10" 3 . Additional Blandford- 
Znajek power could then be important for the spectra of 
ADAFs, both for supermassive black holes in the appar- 
ently quiescent cores of elliptical galaxies (Allen, Di Mat- 
teo & Fabian 1999), and for stellar mass black holes that 
are similarly starved of fuel. 



4 



Blandford-Znajek emission in low luminosity black hole systems 



For the extreme case of isolated stellar mass black holes 
accreting from the local ISM, these results imply that the 
detectability of any rapidly rotating holes is substantially 
better than estimates based on pure ADAF models (Fujita 
et al. 1998). The inferred luminosities are similar to those 
suggested for some unidentified EGRET sources (Dermer 
1997), if we assume (very optimistically) that most of the 
power from the hole is ultimately radiated as gamma-ray 
emission. With better signal to noise observations, vari- 
ability is likely to provide additional constraints on mod- 
els for these EGRET sources. Isolated stellar mass black 
holes, like those in mass transfer binaries, will display vari- 
ability on all timescales exceeding the shortest dynamical 
timescales (of the order of milliseconds) available in the 
system. For holes accreting at high Mach numbers there 



may also be additional quasi-periodic variations arising 
from 'flip-flop' type instabilities in Bondi-Hoyle accretion, 
though the persistence of these features in three dimen- 
sional geometry (Ruffcrt 1997) as opposed to axisymmet- 
ric calculations (Benensohn, Lamb & Taam 1997), is not 
certain. These instabilities have characteristic timescales 
of ~ Ra/coo, the sound crossing time at the Bondi accre- 
tion radius. These variability signatures apply both to the 
high energy emission, and to emission in the optical wave- 
bands that may be detectable via current wide area sky 
surveys (Heckler & Kolb 1996). 

We thank Brad Hansen, Norm Murray, Gordon Ogilvic 
and the referee for useful discussions on these topics. 
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